Engineering fast and stable splitting of matter waves 
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When attempting to split coherent cold atom clouds or a Bose-Einstein condensate (BEC) by 
bifurcation of the trap into a double well, slow adiabatic following is unstable with respect to any 
slight asymmetry, and the wave "collapses" to the lower well, whereas a generic fast chopping splits 
the wave but it also excites it. Shortcuts to adiabaticity engineered to speed up the adiabatic process 
through non-adiabatic transients, provide instead quiet and robust fast splitting. The non-linearity 
of the BEC makes the proposed shortcut even more stable. 
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Introduction. — The splitting of a wavefunction is an 
important operation for matter wave interferometry yj- 
It is a peculiar one though, as adiabatic following, 
rather than being robust, is intrinsically unstable with 
respect to a small external potential asymmetry ^ . The 
ground-state wavefunction "collapses" into the slightly 
lower well so that a very slow trap potential bifurcation 
in fact fails to split the wave except for perfectly sym- 
metrical potentials. An arbitrarily fast bifurcation may 
remedy this but at the price of a strong excitation which 
is also undesired. We propose here a way out to these 
problems by using shortcuts to adiabaticity that speed 
up the adiabatic process along a non-adiabatic route. 
The wave splitting via shortcuts avoids the final excita- 
tion and turns out to be signifficantly more stable than 
the adiabatic following with respect to the asymmetric 
perturbation. Specifically we shall use a simple inver- 
sion method: a streamlined version [dlof the fast-forward 
technique of Masuda and Nakamura [7| applied to Gross- 
Pitaievski (GP) or Schrodinger equations. We have pre- 
viously found some obstacles to apply the invariants- 
based method (at least using quadratic-in momentum 
invariants 0) and the transitionless-driving algorythm 
Q (because of difficulties to implement in practice the 
counter-diabatic terms). 

Fast-forward approach. — The fast-forward method 
0, 0, Q may be used to generate external potentials to 
drive the matter wave from the initial single well to a 
final symmetric double well. The starting point of the 
streamlined version in is the 3-D time-dependent GP 
equation 



(1) 



where the Hamihonian H{t) ^ T + Git) + V{t) in- 
cludes the kinetic energy T, the external potential V , 
and the mean field potential G. Assuming that V is lo- 
cal, (x|y(i)|x') = y(x,t)5(x — x'), it may be written 



from Eq. ([T]) as 

_ ih{^\dMt)) ~ {^\T\m) ~ (xlG(^)IV'(O) 
(x|^(0) 

with (x|7/'(i)) — V'(x, t), whereas 
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VXx,t), 



(x|G(i)|^(i)) = .giV|^(x,t)|2v.(x,t), 

where g is the coupling constant of the BEC and N is 
the number of atoms. For the numerical examples we 
consider ^''Rb atoms, m = 1.44 x lO^^'^ kg. Using in Eq. 
(12) the ansatz 

{^\m) = r(x,t)e'^('''*), r(x,i),<^(x,t) G M, (3) 

the real and imaginary parts of V are 
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where the dot means time derivative. We shall impose 
Im[F(x, t)] = to design a real potential. In addition we 
shall require that the ground state of the initial Hamil- 
tonian H{0) evolves in a time tf into the corresponding 
ground state of the final H{tf), assuming that the Hamil- 
tonian is known at the boundary times. 

In the inversion protocol, r{x,t) is designed first, and 
we solve for in Eq. ([5]) to get Vff '■= Re[F(x, i)] from 
Eq. To ensure that the initial and final states are 

eigenstates of the stationary GP equation we impose r = 
at t = and tf. Then Eq. ([5]) has solutions (t){K,t) 
independent of x at the boundary times Using this 
in Eq. (jH) at t = 0, and multiplying by 6*"^^°^ we get 



— V2+F(x,0)+5|V(x,0)|^ 
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FIG. 1: (Color online) Contour plot of the fast forward po- 
tential Vff in units fku from Eq. Q for the interpolations 
made by (a) Eq. 0, and (b) Eq. Both interpolations 

produce the same initial and final states but b) produces a 
simpler K-shape. Parameter values: uj — 780 rad/s, a = 4 
fim and, tf — 320 ms. 



The initial state 0) is an eigenstate of the stationary 
GP equation with chemical potential — ft0(O) = /-t(O). A 
similar result is found at tf. 

To illustrate this method we consider first a ID lin- 
ear Schrodinger equation (g — 0) and apply the fast- 
forward approach to split an initial single Gaussian state 
r{x, 0) = e"'" ^ /2 (r = rauj / h) into a final double 
Gaussian r{xj,f) = e-r'(^-a)V2 _,. ^-T\x+af/2^ 
vious works [y, [7| use has been made of the interpolation 



r{x,t) = z{t)\ [1 - n{t)]r{x,0) +n{t)r{x,tf) 
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where TZ{t) is some smooth, monotonously increasing 
function from to 1 obeying 7?, = so that r = at the 
boundary times t = and tf, and z(t) is a normaliza- 
tion function. This produces three wave-function bumps 
at intermediate times and a corresponding three- well po- 
tential. Here we use instead the two-bump form 

rix, t) = 2(i)[e-r'(--o(*))V2 + g-r^(.+xo(t))V2]^ (g) 

which generates simpler F-shaped potentials, see Fig. [T] 
We also impose that io(0) = xo{tf) = so f = at 
the boundary times. In the numerical examples we im- 
pose for the Gaussian trajectory the polynomial xo{s) = 
a(3s^ — 2s^), where s = t/tf, and solve Eq. ([5]) with the 
initial conditions = 0) = ^|a;=o = that fix the zero 
energy point. 

Effect of the perturbation. — Now let us assume that a 
small asymmetry affects the splitting process. We model 
this with a potential V\ = Vff + X0{x), where 9 is the 
step function. The splitting becomes unstable, as we 
shall see, but the instability does not depend strongly 
on this particular form, which is chosen for simplicity. 
It would also be found for a linear-in-a; perturbation, a 
smoothed step, slightly different frequencies for the fi- 
nal right and left traps, or a displacement of the central 
barrier 

To analyze the effects of the perturbation A we com- 
pute several "fidelities": The black short-dashed line 
of Fig. [5] represents the structural fidelity Fs = 
lii'o {'^ f)\''Px f))\- It is the modulus of the overlap be- 



tween the (perfectly split) ground state i^oitf) of the un- 
perturbed potential VFF(tf) and the final ground state 
ip^{tf) of the actual, perturbed potential V\. This would 
be the fidelity found with the desired split state if the 
process were adiabatic. Fs{X) decays extremely rapidly 
from 1 at A = to 1/ -^2, which corresponds to the col- 
lapse of the ground state of the perturbed potential Vx 
into the deeper well. 

= lii^o itf)\ip{tf))\, the blue long-dashed line 
in Fig. [21 is the modulus of the overlap between the 
state dynamically evolved with the perturbed potential 
Vx, ^{x,tf) = (x|e^^-*//''|V'(0)), and Vo"(i/), the fi- 
nal ground state of the unperturbed potential VFF{tf ). 
-0(0) = tpx{0) is the initial ground state with Vx{0), but 
the difference with using instead ?/'(0) = i^oi^) in the 
examples shown is negligible, as shown by the overlap 
= li'^x (^)\'^o ~ 1' the green dotted line in 
Fig. 2. 

The flatness of F^"' (A) at small A is in sharp contrast to 
the rapid decay of Fs(X). In practice this feature enables 
us to perform robustly the desired splitting. Note that 
shorter process times tf make the splitting more stable, 
compare the Figs. El^a), (b), and (c). 

Finally, we also calculate Fd = \{i^{t f )\ipx i'^ f))\^ the 
fidelity between the evolved state ip{tf) and the final 
ground state ip^{tf) for the perturbed potential Vx (red 
dotted line of Fig. [2]). For very small perturbations, 
Fo ~ Fs. In this regime the dynamical wave function 
Tpitf) is not affected by the perturbation and becomes 
'ipQ{tf), up to a phase factor, as confirmed also by the 



fact that F, 



(0) 



D 



1 there. We shall understand and quan- 
tify this important regime below as a sudden process in 
a moving-frame interaction picture. As the perturbation 
A increases, the energy levels of the ground and excited 
states of Vx separate and the process becomes progres- 
sively less sudden and more adiabatic. In Fig. [2fc) for 
tf — 320 ms and for large values of A, Fd approaches 1 
again, the final evolved state collapses to one side, and 
becomes the ground state of Vx. For the shorter final 
times in Fig. [UJa) and (b), larger A are needed to make 
Fo approach 1 adiabatically. 

Moving two-mode model. — Static two-mode models 
have been previously used to analyze splitting processes 
Here we add the separation motion of left and 
right basis functions to provide analytical estimates and 
insight. In terms of a (moving) orthogonal bare basis 
0\ fl 
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nian model is 



|i?(t)) 







our two-mode Hamilto- 



H{t) 



A ~6{t) 
-S{t) -A 



(9) 



where S{t) is the tunneling rate [13, Hi and A the energy 
difference between the depths of the two wells [ll]. We 
may simply consider A constant through a given splitting 
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FIG. 2: (Color online) DifTerent fidelities versus the pertur- 
bation parameter A for tlie fast-forward approacli (lines) and 
the two-mode model (symbols). F^'' : (blue) long-dashed line 
and circles; Fd- (red) solid line and squares; Fs- (black) 
short-dashed line and triangles; Fr. (green) dotted line and 
rombs. The (yellow) vertical line is at 0.2/tf. (a) tf = 20 ms. 
(b) tf — 90 ms. (c) tf — 320 ms. The other parameters are 
the same as in Fig. [1] 

process for the time being, and equal to the perturbative 
parameter that defines the asymmetry. A more detailed 
approach that we shall describe later will not produce 
any significant difference. The instantaneous eigenvalues 
are 

E^it)^±^^y + SHt), (10) 

and the normalized eigenstates 

|V^+(t))=sin(f)|L(t))-cos(f)|i?(i)), 

(11) 

l^-(t)) -cos(f)|L(i))+sin(f)|i?(t)), 

where the mixing angle a ~ a{t) is given by tana — 
S{t)/X. 

The bare basis states {\L{t)), \R{t))} are symmetrical 
and orthogonal moving left and right states. Initially 



FIG. 3: (Color online) Coordinate representation of the time 
dependent bare basis. The right vector \R{t)) is plotted for 
the parameters of Fig. [1] The left vector \L{t)} satisfies 
{x\R{t)} = {—x\L{t)) and is orthogonal to it. The side (neg- 
ative) peak eventually dissapears. 




FIG. 4: (Color online) The same fidelities as in Fig. [2]versus 
the perturbation parameter A using the fast-forward approach 
for a Bose-Einstein condensate. F^^ : Blue long-dashed line; 
Fd: red line; Fs: black short-dashed line. Parameter values: 
tf — 320 ms, gN / {TiLijaho) ~ 1.38. The rest are the same as 
in Fig. [H {aho = y^h/{muj)). 

when they are close enough (and (5(0) >> A), the in- 
stantaneous eigenstates of H are close to the symmetric 
ground state \ipo i^)) = Tf (1^(0)) + 1-^(0))) and the an- 
tisymmetric excited state \^o{0)) = -j^{\L{0)) - \R{0))) 
of the single well. At i/ we should distinguish two ex- 
treme cases: i) For S{tf) >> A the final eigenstates of 
H tend to = ^mtf)) ± mtf))) which cor- 

respond to the symmetric and antisymmetric splitting 
states, ii) For S{tf) << A the final eigenfunctions of 
H collapse and become right and left localized states: 
|V'^(t/)) - \L{tf)) and l^+itf)) = \R{tf)). Since 6(1 f) 
is set as a small number to avoid tunnelling in the final 
configuration, the transition from one to the other regime 
explains the sharp drop of Fs at small A « S{tf). 

Dynamics of the two-mode model. — We define a 
moving-frame interaction-picture wave function — 
A^^, where A = E,3=L,i?J^(*))(/3(0)| and is the 
Schrodinger-picture wave function. In principle it obeys 
ihijj^ = {Ha - Ka)iP^, with Ha = A^HA, and Ka = 
ihA^A but, for real {x\R{t)) and {x\L{t)), the symmetry 
{x\R{t)) = {-x\L{t)) makes Ka = 0. 

We may invert Eq. ((TT|) to write the bare states in 
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terms of the ground and excited states and energies, and 
get S{t) from the energy sphtting in Eq. ([10]). The 
two-level model approximates the actual dynamics by 
first identifying and E^{t) with the instanta- 

neous ground and excited states and energies of the un- 
perturbed fast-forward Hamiltonian. We combine them 
to compute the bare basis in coordinate representation. 



/3'/3 
A ' 



Then we compute the matrix elements (/?' | H\ \/3) = H 
{Hx = T + Vx) for 13^13'. By comparison with 
Eq. © we get S{t) ^ -2H^^/h = -2H^^/h. For 
the diagonal, j3 = /?', wc may for consistency calcu- 
late X'{t) := 2(iJf« - Va)/h ^ -2(i?f ^ - Vo)/h, where 
Vo = Vo{t) = [E^{t) + E+{t)]/2 is a shift to match the 
zero energy point between the fast-forward and the two- 
mode model. A' differs slightly from the constant A at 
short times. In our numerical calculations the results 
of substituting A by A' are hardly distinguishable in the 
final fidelities so the treatment with A is preferred for 
simplicity. Adjusting from the fast-forward approach the 
values of 5{t) we solve the dynamics in the moving frame 
for the two mode Hamiltonian. The initial state may be 
the ground state of the perturbed or unperturbed initial 
potential, as for the fast-forward calculations the results 
can hardly be noticed in the figures. The comparison 
with the exact results are excellent, see the symbols of 
Fig. [21 so the two-level model provides a powerful inter- 
pretative and control tool. For further insight we shall 
perform further approximations. 

Sudden and adiabatic approximations. — The fidelities 
at low A may be understood with the sudden approxima- 
tion flS*!. Its validity requires 



tf < H/AHa, 



(12) 



where AHa = WO)\Ha |^A(0)) - (V'(0)|i?A|V'(0))2. 
We take |V'(0)) = |V'o"(0)), and i?! = ^ J^' dt'HA{t'), 
where the matrix elements of HA{t') in the bare basis 
{1/3(0))} coincide with the matrix elements of H in Eq. 
©, when the later are expressed in the basis {\P{t'))}- 
The condition for the sudden approximation to hold be- 



(13) 



A vertical line marks 0.2 /tf in Fig. 2. 

The subsequent increase of F]j for increasing A can 
be explained using the complementary adiabatic approx- 
imation. The adiabaticity condition is here IJ, Il5| 



Ki^ximi^m « jiiExit) - Et{t)\. (14) 

Taking into account Eqs. pH)) and pT|) it takes the form 



XS{t) 



2(A2+5(i)2)3/2 
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FIG. 5: (Color online) Fidelities for the Bose-Einstein con- 
densate, same line codes as in Figs. 2 and 4. \/(hw) = 0.02, 
a = 4 /im, tf = 320 ms. 



Thus the structural and adiabaticity criteria delimit the 
small and large A values for which Fd ~ 1. 

Bose Einstein condensates. — We use now the stream- 
lined fast-forward method to analyze the effects of the 
asymmetric perturbation in the splitting of a BEC. There 
is no analytical expression for the ground state of a BEC 
in a harmonic trap, so to mimic Eq. ([8]) we get first nu- 
merically the ground states xn{x) and Xii-i^) of a ti^r- 
monic trap that holds a BEC with TV and N/2 particles 
and design the time evolution that connects these two 
states by interpolation as f{x,t) = [1 — TZ{t)]xN{x) + 
n{t)xiL{x), where 7^(^) = Ht/tff - 2{t/tff. Finally 
r{x, t) is constructed by displacing and summing these 
functions, 

r{x,t) = -l-^\^f[x - xo(t),t] + J[x + xo{t),t]^, (16) 

where z{t) is a normalization factor and xa{t) = aTZ(t). 
Note that this form reproduces Eq. ([S]) for iV = 0. We 
then get Vff from Eqs. ([4]) and ([5]) and evolve in a 
perturbed potential Vx{t) = Yppit) + XO{x) the initial 
ground state of the CP equation associated with Va(0). 
The same fidelities as in the linear case may be computed 
and the results are shown in Fig[31 The structural fidelity 
(black short dashed line) is not much different from the 
linear case, i.e., splitting by adiabatic following is also 
very unstable for the condensate. Similarly, the dynam- 
ical fidelity Ed (red solid line) drops abruptly at small 
A in a sudden regime to increase more slowly later on 
towards the adiabatic regime. There is however a very 

(0) 



< 1. 



(15) 



remarkable stabilization of F^ ' « 1 (blue long dashed 
line) with respect to the linear dynamics, as seen in Fig. 
[3] with respect to A (the scale for the A-axis differs by 
an order of magnitude in Figs [5] and [1]) and Fig. [S] with 
respect to the non-linear coupling constant. In Fig. [5| 
we see that the relevant fidelity saturates to one. This 
means in summary that fast splitting via a designed fast- 
forward potential is significantly more robust versus the 
asymmetry for the condensate, as the non-linear term 
compensates for the external potential asymmetry. 
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Discussion. — We have designed simple F-shaped (po- 
sition and time dependent) potential trap bifurcations to 
split matter waves rapidly without final excitation, avoid- 
ing the intrinsic instability of the adiabatic approach with 
respect to slight asymmetries. Incidentally, we also avoid 
or mitigate in this manner the decoherence effects that af- 
fect slow adiabatic following. The bifurcation may be ex- 
perimentally implemented by means of spatial light mod- 
ulators A simpler approximate approach would in- 
volve the combination of Gaussian beams. Further stan- 
dard manipulations may be combined with the proposed 
technique, in particular a differential phase among the 
two final parts may be imprinted by illuminating one of 
them with a detuned laser. 
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